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We present an analysis of the mass spectrum of heavy quarkonia with non-degenerate quark masses. The heavier 
(bottom) valence quark is treated in a non-relativistic fashion and the other (charm) is a relativistic Wilson-like 
quark using the improved SW action. Such states provide and interesting probe between the relativistic B meson 
states and the non-relativistic bottomonium states. 



1. Introduction 

The mass spectrum of quarkonia with valence 
bottom and charm quarks are interesting in a 
number of aspects. Though only a handful of 
events have been detected, upcoming B-factories 
should produce copious amounts of them. Hence, 
we are put in the interesting circumstance of mak- 
ing a pre-diction rather than a post-diction for 
non-exotic hadron masses. As they interpolate 
the non-relativistic bottomonium states and fully 
relativistic heavy-light states such as the B and 
Bg they provide an interesting probe of QCD dy- 
namics. Of course, the calculation of these masses 
are also a first step in the calculation of decay con- 
stants and semi-leptonic widths, and these are of 
interest phenomenologically as an accurate cal- 
culation of /bc is crucial to obtaining as accurate 
as possible a result for fs from LEP data j|] and 
the semi-leptonic decay Be J/ ipW will give an- 
other estimate of the matrix element | Vch \ , which 
will be seen at the LHC-B g]. 

2. Lattice method 

Previous quenched lattice calculations of the 
bc spectrum has been carried out at /3 = 
5.7, where both valence quarks are simulated us- 
ing the non-relativistic formalism On such 
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coarse lattices, arrib and arric > 1 and the effec- 
tive Green's functions for both quark masses can 
be calculated without introducing large radiative 
corrections to the hamiltonian. However, it has 
been demonstrated that for such lattice spacings 
the charmonium spectrum is not reproduced very 
accurately. For quenched gluon fields at /3 = 6.2, 
quenched one can still use the non-relativistic ap- 
proach for the b quarks and the relativistic for- 
malism can now be used for treating the c quark. 
As one expects the dominant errors to be a func- 
tion of arric, which still approaches 1, we use the 
Sheikholeslami-Wohlert action , using a mean- 
field estimate of the clover coefficient Cgw • While 
not eliminating all errors up to ©(a^), this mean 
field estimate approach does reduce errors signif- 
icantly g. 

3. Computational Details 

On a 24^ x 48 lattice at /3 = 6.2, correspond- 
ing to a lattice spacing of approximately 2.7 GeV, 
two point functions were constructed for the fol- 
lowing states; (l,2)i5'o, (1,2)351, and3Po,i,2 
(using the notation (i?)^''+^Lj, where R is the 
radial excitation). The operators were calculated 
using Coulomb gauge fixed smeared sourcesQ. 
In order to maximise the statistics, the non- 
relativistic Green's functions were calculated for- 
wards and backwards in time, so that as many 
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Table 1 

Mass parameters used for states 



K 




# Configs. 


0.132 


1.22 


100 


0.126 


1.22 


100 


0.126 


2.44 


100 



of the componants of the relativistic propagator 
could be used. In table |l| the quark mass param- 
eters and the number of configurations are listed. 
We note that k = 0.126 and arrib — 1.22 corre- 
spond roughly to charm and bottom quark masses 
from studies of charmonium and bottomonium 
systems respectively. The NRQCD hamiltonian 
was corrected to 0{mv'^)^ using tadpole improved 
estimates of the relevant coefficients. As the 
configurations are quenched, the lattice spacing 
should be determined from one of the be observ- 
ables. As the splitting Smp^s = ^ Pi — ^S' for char- 
monium and bottomonium only changes by a few 
MeV, we took an average of these splittings and 
hence assumed that Smp^sipc) = 0.4545 GeV. 

3.1. Mixing 

One expects that the ^Pi and '^Pi states to mix 
(this has been observed at /? = 5.7 Q and with 
two relativistic propagators at higher /3's with dif- 
ferent masses around that of charm). Therefore 
we calculated the following matrix of correlators 

< iPi^P/ > < iPi^p/ > 
<3Piip/> <3Pi3p/> 

However, for all three combinations of quark 
masses, the off-diagonal elements of this matrix 
were very noisy and were consistent with zero. 

4. Results 

In figure |] the results for the be spectrum us- 
ing the the best estimates for the bottom and 
charm quark masses are compared with the pre- 
vious lattice calculation at /3 = 5.7 (|] and a po- 
tential model calculation The R = 1 states 
at both /3's are consistent with each other. For 
the S'-wave states, where we have a measurement 
of the R = 2 states, the agreement is to within 
two standard deviations, albeit with large errors. 



The potential model calculation agrees with the 
lattice results except for the hyperfine splitting. 
It is not clear whether this difference is due any 
systematic error in the potential model or by a 
choice of the scale using Smp^s- 
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Figure 1. A comparison of the results at /? = 
5.7[§, /3 = 6.2 [k = 0.126, arub = 1.22) and a 
potential model calculation. The labels on the 
horizontal axis are arbitrary. 



In figure || the results for the be spectrum for 
different combinations of quark masses at /? = 6.2 
are shown. Once again, there is a general broad 
agreement. For the splittings, we find that the 
splittings behave as expected. For example, in fig- 
ure H we see that as the b quark mass is increased, 
and the system behaves more like a heavy-light 
bd state the splitting tends to zero. Likewise, the 
fine splitting (in figure |^) is constant, indicating 
that Smp^s should be insensitive to variations in 
the quark masses. 
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Figure 2. A comparison of the results at /3 = 6.2 
for the different combinations of quark mass. The 
labels on the horizontal axis are arbitrary. 



Figure 3. A comparison of the ^Si — ^Sq split- 
ting at (3 — 6.2 for the different combinations of 
quark mass. The labels on the horizontal axis are 
arbitrary. 



5. Conclusions 



6. Acknowledgements 



It is very promising that the quenched lattice 
results at two different lattice spacings are in good 
agreement with each other, despite potential con- 
cerns over the size of the discretisation error at 
(3 — 5.7. The broad agreement with potential 
models is also quite interesting as it is not clear 
that a semi-relativistic system such as this could 
be modelled with a simple potential. Given the 
amount of research that has been done on the bb 
and b{u, d) spectrum, it would be interesting to 
see how the splittings of these states vary as the 
masses of the quarks are changed. 

The absence of mixing between the "^Pi and ^Pi 
states at f3 — 6.2 is quite strange, given its pres- 
ence a.t /3 — 5.7. It is of course possible that the 
mixing is a strong function of the lattice spacing 
and that the mixing is very small in the contin- 
uum. This would give a qualitatively different 
result from potential models. In order to verify 
if this is true, the calculation should be repeated 
using some other set of operators, defined in a 
different gauge or a gauge invariant fashion. 



The numerical work for this presentation was 
carried out on a Cray J90 at the EPCC in Ed- 
inburgh. HPS is supported by a grant from the 
Leverhulme foundation. 
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Figure 4. A comparison of the ^Pq — ^Sq split- 
ting at f3 = 6.2 for the different combinations of 
quark mass. The labels on the horizontal axis are 
arbitrary. 



